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Abstract 

A mean-field theory based on Gibbs-Bogoliubov inequality is constructed to study 
the interactions between two like-charged polyions. It is shown that contrary to 
the previously established paradigm, a properly constructed mean-field theory can 
quantitatively account for the attractive interactions between two like-charged rods. 



One of the most fascinating problems that has recently appeared in the field of 
condensed matter physics is the discovery of attraction between like-charged 
macromolecules [1]. This attraction plays a fundamental role in various bio- 
logical processes such as the condensation of DNA [2,3] and the formation of 
fibers composing cellular cytosceleton [4] . The attraction between like-charged 
colloids has also been observed in various experiments and simulations [5-8] . It 
has been noted that the attraction appears only in the presence of multivalent 
counterions. 

A number of models have been proposed to try to explain the mechanism 
of these strange phenomena. It is now clear, from both simulations and ex- 
periments, that this effect is purely electrostatic and is produced by strong 
many body interactions present in polyelectrolyte solutions. In a beautiful set 
of experiments Tang et al[9] demonstrated how addition of simple monova- 
lent salt produced dissociation of the actin bundles. The F-actin chains are 
highly charged polymers, which inspite of their large negative charge density, 
aggregate in well defined bundles in the presence of polyamines. However, this 
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bundling can be reversed by addition of simple monovalent salt which screens 
the electrostatic interactions between the polyions and the multivalent coun- 
terions. 



The first explanation of attraction between like charged surfaces in the pres- 
ence of multivalent counterions was advanced by Kjellander and Marcelja [10] 
based on the integral equation formalism. From the numerical solution of the 
AHNC equation these authors came to conclude that for sufficiently high sur- 
face charge, an attraction can arise between like charged plates. A very sim- 
ple physical picture to explain the mechanism of attraction was advanced by 
Rouzina and Bloomfield [11], and extended by Shklovskii [12]. These authors 
proposed that the condensed counterions around the two plates form strongly 
coupled Wigner crystals. In the case of rod-like polyions, a similar explana- 
tion has been advanced by Arenzon et al on the basis of an exactly solvable 
model[13,14]. A different mechanism, relying on correlated fluctuations, has 
been proposed by Ha and Liu [15], but has been criticized by Levin et al [16]. 

Since the beginning of the study of this interesting phenomenon there has been 
a general consensus that the attraction must arise as a result of correlations 
of condensed counterions [17]. It was, therefore, implicitly assumed that no 
mean-field theory would be able to account for this phnenomenon. This belief 
was further reinforced by the solutions of Poisson-Boltzmann equation (PB) 
which, of course, did not predict any attraction. Not all mean-fields, however, 
are equal. In this paper we shall present a mean-field theory, which quantita- 
tively accounts for the attraction between like-charged rods in the presence of 
condensed multivalent counterion. 



We consider two parallel polyions modeled as rigid rods, each having Z charges 
of value —q, spaced uniformly with separation b along the length. The rods 
are separated by distance d = xb. The strong electrostatic interaction between 
the polyions and the multivalent counterions present in solution leads to coun- 
terion condensation [18-21]. The effect of n, a-valent condensed counterions, 
is approximated by the renormalization of local charge. Thus, if one of the 
charged sites of a polyion has an associated condensed counterion its effective 
charge becomes —q(l — a). Note that in this simple model the condensed coun- 
terions are assumed to reside only on top of the charged sites. The net charge 
of each polyion is (Z — an)q. The Hamiltonian for the interactions between 
the two rods is [13], 

2^^^ r(i,m;i',m') 



where we have introduced to — 0, 1 to label the two polyions. The distance 
between two charged sites of the polyions, (i, m) ^ (i', m') , is r(i, to; i', to') = 



2 



by\i — i'\ 2 + (1 — 5 mm ')x 2 , and a™ is an occupation variable such that a™ = 1 
if the i'th site of m'th polyion has an associated counterion and cr™ = if this 
site is unoccupied. 

The mean field theory can be constructed with the help of the Gibbs-Bogoliubov 
bound for the free energy [21], F < T = F + < 7i — 7i Q > . The average, 
< . . . >o, is performed with respect to the trial Hamiltonian 7{q. To make the 
calculation as simple as possible we shall take this to be of one body form, 

n = - q ^-^r)vT, (2) 

i,m 

where ^ is the mean electrostatic potential experienced by the i-th monomer 
of the rod m. The upper bound for the free energy can now be calculated, 



jr = J_Y Y g 2 (l-<)(l-<') 
' 2^,fc lm ^ =0 r(i,m;i>,m>) 

+ kE [(1 - a<) ln(l In <] , (3) 

P i m=0,l 

where the average occupation per site is n m = (o\ m ) and the constraint, 
J2i n T = n i 1S implicit. The optimum upper bound is obtained from the mini- 
mization of the functional in Eq. (3). We find, 

™ n 

(4) 



n + £,-(! - nf) expfa^^ - (j>f)] ' 



where ^ = q 2 /DksTb and 



„ an™ - 1 „ an) m - 1 

C = £Thr + £ I ' ■ (5) 



Here m is the reduced electrostatic potential experienced by the condensed 
counterion of the i'th site of the m'th polyion. These equations can be solved 
numerically, producing the positional distribution of the condensed counteri- 
ons on the two polyions, Fig 1. We make the fundamental observation that the 
two profiles are not equal. Thus, the mean-field theory breaks the symmetry 
between two polyions! This is clearly an artifact of mean-field approximation. 
Obviously if the density profiles would be calculated exactly by an explicit 
solution of the partition function, they would be identical. There is no way of 
breaking the symmetry between two identical finite sized polyions. In the case 
of exact solution, there would, however, exist very strong correlations between 
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the condensed counterions on the two polyions. These would provide an im- 
portant contribution to the total free energy. Since the mean-field theory does 
not account for these correlations, in order to establish an optimum bound, it 
breaks the symmetry between the two rods. 

The horizontal component of the force between the two polyions is, 

q 2 v (l-an?)(l-anj) x 
k DPjj x* + \i-j\* v / x2 + |,_ J | 2 ' W 

where i and j correspond to the sites on rods and 1, respectively. The density 
profiles n° and n 1 are obtained from the solution of Eqns. (4,5). For short 
separations between the polyions the force becomes attractive. This is the 
result of the symmetry breaking discussed above. The force calculated using 
the mean-field theory, Fig. 2, is in quantitative agreement with the Monte 
Carlo simulations and the exact solutions. 
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Fig. 1. Density profile for two rods (solid and dashed lines) for Z = 20 and n = 9 
divalent counterions, x = 0.4. Notice the staggered configurations along the two 
polyions. 

The attractive force is short ranged and appears only if the number of (a > 2) 
counterions is larger than a threshold, n = Z/2a. For a — 1 the force is always 
repulsive, which is in full agreement with the experimental evidence on the 
absence of attraction if only monovalent counterions are present [2]. 
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Fig. 2. The horizontal component of the force between two like-charged rods. Z = 20 
and n = 9, the net charge on each rod is —2. The points are obtained using Monte 
Carlo [13], while the solid line is the mean field result. 
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